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We analyse the use of entangled states to perform quantum computations non locally among 
distant nodes in a quantum network. The complexity associated with the generation of multiparticle 
entangled states is quantified in terms of the concept of global cost. This parameter allows us to 
compare the use of physical resources in different schemes. We show that for ideal channels and 
for a sufficiently large number of nodes, the use of maximally entangled states is advantageous over 
■ uncorrelated ones. For noisy channels, one has to use entanglement purification procedures in order 

' to create entangled states of high fidelity. We show that under certain circumstances a quantum 

network supplied with a maximally entangled input still yields a smaller global cost, provided that n 
belongs to a given interval n £ [n m i n ,n max \. The values of n m i„ and n max crucially depend on the 
purification protocols used to establish the n— processor entangled states, as well as on the presence 
of decoherence processes during the computation. The phase estimation problem has been used to 
illustrate this fact. 



I. INTRODUCTION 



Consider a quantum computation which can be divided into subroutines so that each subroutine can be run on a 
separate quantum processor. The processors may be placed at different locations/nodes of a computational network 
and may exchange data with a selected central processor (!]] . Each processor operates on a partial input which has a 
fixed size. The partial inputs may be independent of each other, correlated or even entangled. When the computation 
is finished, the central processor, after collecting partial outputs from the other processors, stores the global output. 
This type of distributed computation may be repeated several times to yield a desired result and as such it features 
frequently in quantum parameter estimation procedures e.g. the phase estimation in frequency standards (2|^]. In 
some computational tasks, e.g. estimating a given parameter with a prescribed precision, the number of repetitions 
depends on the form of the input state — some entangled states require less repetitions than uncorrelated inputs. In 
the case of correlated input states, we have to pre-compute the input state for each run of the computation and this 
involves an additional use of physical resources. Are we still better off when the complexity of the pre-computation 
• t-h is included? How shall we include and compare the use of different physical resources? 

In this paper we quantify this complexity by introducing the notion of cost of physical operations, such as the cost 
of establishing an entangled pair over a channel, the cost of transmitting one classical bit between components, the 
cost of running a quantum processor, etc. and discuss the performance of the distributed quantum computation when 
the inter-processor quantum communication is prone to errors, i.e. when the quantum channels among the processors 
are noisy. 

This paper has been organized as follows. In section II we introduce the notion of cost in distributed quantum 
computation and show that if the ratio between the number of repetitions for the entangled and separate inputs 
decreases fast enough with the size of the network then, above some critical size, the computation can be made 
cheaper using entangled inputs. This behaviour is illustrated explicitly in the case of estimation of a small phase shift 
for qubits, which was considered in Ref. jjj. In section III we introduce the phase estimation problem for disentangled 
and maximally entangled states and show how to compare the two scenarios. In section IV we analyse the ideal case of 
noiseless channel and error-free computations. Noise along the channels linking the nodes of the network is taken into 
account in section V and the use of different purification schemes is considered. We show that for certain purification 
protocols it is advantageous to use entangled states, whereas for others the cost of the precomputation is not off-set by 
the subsequent reduction in the number of repetitions. In section VI we analyse the effects of decoherence during the 
computational process at each node and show how they affect the results obtained in the preceding sections. Finally, 
in section VII we summarise the main results of this work. 
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II. COSTS OF DISTRIBUTED QUANTUM COMPUTATION 



Let us start out by considering a generic scenario for distributed quantum computation: a central processor A and 
7i — l processors labelled Bj (i = 1, 2, . . . , n — 1) represent the nodes of a quantum network. These nodes agree on 
performing a given computation which consists of three steps: (1) Precomputation: In order to prepare the initial 
state of all the nodes they exchange certain classical and quantum information; (2) Computation at each node: Each 
of the nodes performs a well defined operation locally, followed by a measurement; (3) Communication of the results: 
The B nodes report the outcomes of their measurements to the central node by sending k bits of classical information. 
With this information, the central node estimates the outcome of the complete computation. The computation gives 
the correct result with certain probability, and therefore it has to be repeated a number of times in order to achieve 
a prescribed precision. 

We are interested in the cost of the computation in terms of the number of uses of the processors, and the amount 
of classical and quantum communication involved in the whole computation. Let us denote by P(n) the cost of the 
precomputation, i.e., the cost of establishing the initial state for the n nodes. In general, if one wants to create 
an entangled state of all the nodes one has to send quantum information through the channels. Besides, due to 
the presence of noise during transmission and processing, the nodes will have to use either purification or error 
correction methods, which will require in addition some classical communication. Thus, the cost of precomputation 
will depend on the costs of sending qubits and classical bits through the communication channels that link the nodes. 
We denote by Z the cost of running a quantum processor at each node, and by Y the cost of sending the outcomes 
of the measurement from one node to the central node. Finally, we will denote by R(n) the number of times the 
computation has to be repeated in order to obtain a prescribed precision. With these definitions, one can calculate the 
total cost of the computation. We wish to analyze the advantages of using entangled input states in the computation 
with respect to the case of initial uncorrelated states. Thus, we consider the following two scenarios: 

(1) Disentangled states: If the initial state of the processors is disentangled, no communication is required in the 
precomputation. Therefore we take P(n) = B and obtain 

Ci(tj) = Ri(n)[nZ + (n- 1)Y]. (1) 

(2) Entangled states: If the initial state of the processors is entangled, communication is required in the precompu- 
tation. We have 

C 2 (n) = R 2 (n)[P 2 (n) +nZ + (n - 1)Y}. (2) 

We can now evaluate the ratio between the cost C 2 (n) corresponding to entangled inputs and the cost Ci(n) for 
independent processors. We obtain 

C 2 {n) = R 2 (n) P 2 {n) + (n - 1)Y + nZ 

Ci(ra) _ Ri{n) (n-l)Y + nZ ' ^ ' 

This ratio depends crucially on the ratio between the repetitions needed in each case, as well as on the cost of the 
precomputation. Under ideal conditions, the use of entangled states in general decreases the number of repetitions 
required, i.e. R 2 (n) < Ri(n). The use of entangled states will be cost efficient if what one gains in the number of 
repetitions compensates what one loses in the precomputation. Thus, we expect that there will be a certain n min 
such that if n > n m i„, the use of entangled states is cost efficient. On the other hand, in the non ideal situation 
in which one has noise either in the quantum channel or during the computation, R 2 (n) may increase with n more 
rapidly than R\{n) since entangled states are more prone to errors. Thus, we would expect that for a specific task 
there is a maximum value n max such that for n < n ma x the use of entangled states is cost efficient. Therefore, if 
n m in < n max there will be an interval [n m i n , n max ] in which entangled states are more advantageous. In the opposite 
case n m i n > n max , it will be more convenient to use disentangled states. The values of n m i n and n max will depend on 
the specific computation and of the methods used to create the entangled states. In the following we will concentrate 
on the specific but relevant problem of phase estimation addressed in Ref. |]J . We will first analyze the ideal noiseless 
case. Then we will consider the presence of noise in the quantum channel, and analyze two ways of overcoming this 
noise using entanglement purification. Finally, we will consider the effects of decoherence during the computation. 

III. THE PHASE ESTIMATION PROBLEM 

As an illustrative and, in spite of its simplicity, important example, consider a network of n processors each 
performing computation C defined as a small conditional phase shift on a qubit 
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|0>— |0), (4) 
>e**|l). (5) 

Computation C is performed at each of the n nodes (A, Bi). Each run consists of a conditional phase shift and the 
subsequent measurement protocol. Computation C is then reset after each repetition: we assume that no extra-phase 
accumulation is allowed by means of consecutive runs of the computation C on the same qubit before the measurement 
is performed. 



A. Disentangled states 



Without inter-node entanglement, the best way to estimate 4> is to prepare each node in the initial state 

l*) i =^(|0) i + e- i * l |l) i )- (6) 

where (f>i is a given phase that can be adjusted from computation to computation. Computation C is then applied, 
followed by a Hadamard transformation H given by H = (|0)(0| - |1)(1| + |0)(1| + |1)(0|). The last step is the 
independent measurement of each qubit in the computational basis. The result of the measurement will be either or 
1 with probabilities p% and 1 — p\, respectively. Each of the Bi nodes then transmits one classical bit, corresponding 
to the result of the measurement, to the central node A. This process is repeated R± times yielding a binomial 
probability distribution. In this way one can estimate <j> at node A with precision 



£ i = ,d P1 , m > ( 7 ) 



where Api = pi(l — pi) is the variance of the binomial distribution. In general, e\ will depend on <j> and tf>\. As soon 
as the first outcomes of the measurements are obtained (first repetitions) the value of <f>\ can be adjusted in order to 
minimise t\ H. 



B. Entangled states 



Let us now assume that the initial state of the n nodes is an entangled state of n qubits of the form 

| * 4d ) = -L(| 000 ... 0) + e- m ^ | 111 ... 1)) . (8) 

Under ideal conditions, this state can be obtained as follows. The central processor at the central node generates n— 1 
EPR pairs and sends one member of each pair to the remaining nodes. An EPR pair shared between node A and Bi is 
referred to as the A — Bi pair. In order to obtain the state in Eq. (^|) we pick up one of the n — 1 qubits at node A and 
using it as a control qubit we apply the quantum controlled-NOT operation | ei) | £2} — > |ei)|ei © £2) (ei.2 = 0, 1 and © 
denotes addition modulo 2) with the remaining n — 2 target qubits at node A. Then we measure the n — 2 targets in 
the computational basis. At this stage we have already established an entangled state of all n nodes; in order to put 
it into the form (^) we simply perform operation NOT (a x ) at location Bk if the result of the measurement performed 
at A on the qubit belonging to the A — Bk pair was 1. Finally, all the nodes perform a phase shift transformation 
with angle <j)2- In Figure 1 we have depicted a set-up for the simplest case, involving only three processors. Once 
EPR pairs of the form -^=(\ 0) A \0) B . + | l) A 1 have been established between nodes AB\ and AB 2 (via channels 

represented in the figure by a thick line networking the central node with the other two), the central node A performs 
a CNOT operation between the two qubits stored in A. A measurement of the target bit in the computational basis 
reduces the three-node composite state to either the state -^(| 0) A | 0) Bi | 0) Ba + | 1) A | | l) Ba ) (outcome 0) or to 

the state ^(| 0) A | 0) Bl | + | 1) A \ l) Bi | 0) Bo ) (outcome 1). In the latter case, node B 2 has to invert the state of 
its qubit by means of a local operation. Therefore, the pre-computation requires nodes A and Bi to exchange classical 
information, as illustrated in the figure by the dotted line connecting those nodes. 

Once a state close to the ideal state (||) is prepared, at each node we run computation C followed by the Hadamard 
transform. A measurement on the computational basis is then performed at each node. Nodes Bi report their outcomes 
to node A by broadcasting one bit of information and the overall parity of the reported bits and the outcome at node 
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A is calculated at A. This will give bit value or 1 with probabilities p 2 (<^) an d 1 — P2 (</>)■ The procedure is repeated 
R 2 times and gives estimation of cp with precision 



^P2 

r R2 



V ^P2 ( C\\ 

£ 2 = |dpa , > ( 9 ) 



where Ap 2 — p 2 {^ — P2)- In general, t2 will depend on and cf> 2 . As soon as the first outcomes of the measurements 
are obtained (first repetitions) the value of <j> 2 is adjusted in order to minimise e 2 |5|]. 



C. Comparison 

In order to compare the two procedures, we impose that the precision required is the same, i.e. t\ = 62 = e. We 
obtain 



R 2 (n) _ Ap 2 ( dpx\ 2 ( dp 2 s 
Ri{n) ~ ~A~p[ 



(10) 



where in this expression the phases <pi,2 are assumed to be chosen independently of each other in order to minimize 
e. Once the values of Api and dpi/dcj) are known, one can substitute this expression in (§). The value of P 2 {n) will 
depend on the purification procedures used in the generation of the state (£jh . 

IV. IDEAL CHANNELS AND COMPUTATIONS 

We consider first the simple situation in which no decoherence is present. In this case 

px = -[l-cos(^-^i)], (11a) 

p 2 = ~[1 - cos(n<p - ri(p 2 )]. (lib) 

One obtains e\ = 1/ (nRi) 1 / 2 and e 2 = ^/(nR 2 ) 1 ^ 2 independently of the values of </>i,2- Therefore, we have R\{n) = 
l/(ne 2 ) and R2(n) = l/(ne) 2 . On the other hand, the cost of the precomputation using the procedure described above 
is simply — ( n ~ l)-^ + (n — 2)Y", where X is the cost of sending one qubit from the central node to any other 

node. We finally obtain 



C 2 1 / (2n - 3)Y + (n- l)X + nZ 



d n\ (n- l)Y + nZ 



(12) 



This expression implies that for n larger than a certain value n m i n , the global cost for computation with entangled 
states is smaller than the one with independent states. It can be easily checked that for (2Y + X + Z)/(Y + Z) ^S> 1 

_ 2Y + X + Z 
n min ~ Y + Z 

If the cost Z is much smaller than X and Y, and Y is much smaller than X, the threshold value is given by the ratio 
of the costs of distributing entanglement and classical communication. Figure 2 illustrates this behaviour. 



V. NOISY CHANNELS 



We have seen in the preceding sections that the use of entanglement is cost efficient above a certain threshold in the 
number of nodes of a quantum network. However, this result holds under the assumption that the channels networking 
the nodes are ideal. In particular, this implies that ideal entangled states | ^>id) can be distributed among the n nodes. 
In reality, this will never be the case, and therefore one has to analyze what will happen for noisy channels. While 
creating the state \^id) (using, for instance, the protocol exemplified in Fig. 1), there will be errors. The state 
will no longer be a pure state, but will rather be described by a density operator p ^ \^id){^id\- The closeness of 
this state to the ideal one is measured by the fidelity F = id\p\^ id) ■ This means that the number of repetitions 
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R,2(n) required to perform the computation to a prescribed precision will increase (since the probability distribution 
of obtaining the right outcome becomes worse). On the other hand, one may use entanglement purification in order 
to increase the value of Fq and therefore to reduce the number of required repetitions. However, this requires a 
higher precomputation cost Pzin). In this Section we analyze this problem for two different purification protocols, 
specially suited for different situations. In order to focus on the role of noise along the channels, we will assume that 
all local operations (both the ones required for the establishment of the entangled states and the ones involved in the 
computation at each node) are error free. 

We assume that using a two-qubit purification protocol and after s steps one creates an entangled pair between 
nodes A and Bi with fidelity F s < 1. Let us denote by P\(s) the cost required to create such a state. Once we have 
the entangled pairs we use the method described above (see Fig. 1) to create the entangled state among the n nodes. 

Assuming that the pairs arc in a Werner-like state, the fidelity of the n-qubit state will be F n ~ In order to 

estimate how the results are affected by noise, we consider for simplicity 

p = Xn \^ ld )^ ld \+ l —^I l (14) 

where F n = x n + (1 — x n )/2 n . If we perform the computation with this state instead of the ideal one we obtain 

p 2 = - cos(n4> - nfa)] + — (15) 
Therefore, to estimate the parameter <f> with precision e, the computation has to be run a number of times 

l L i 

x^ sm 2 [ncj) — nfa) 

which now depends on (f> and <f>2- We take the optimal choice of </>2, which gives 



Mn) = Z j ns [l+ _ 9 ,_,_„" J (16) 



Mn) = -o-Ur , (17) 



and obtain the cost ratio 



C 2 (n) _ 1 /(n-l)P (s) + (2n-3)F + nZ 
din) ~ ^xj \ (n-l)Y + nZ 



(18) 



Taking into account that x n ~ F n ~ .Fj" 1 we see that for a sufficiently high value of n, this ratio will be as big 
as we please. This implies the existence of a value n max such that for n > n max there is no gain in using entangled 
states. In fact, it may happen that there is no gain for any value of n. On the other hand, the value n max will 
depend on Pq{s) and F s , i.e. on the specific purification procedure. Besides, as before, there will be a value n m i n such 
that if n £ [n TO , n , n max ] then entangled states can be cost efficient. We will illustrate these features for some specific 
purification protocols. 



A. Purification Scheme 1 



Let us assume that we can create pairs between node A (central node) and node Bi (i = 1, n — 1) with fidelity 
Fq. For simplicity we assume that they are in a Werner state. We consider that local operations are perfect, and 
therefore one can use the purification procedures of Refs. 1^,0). One can easily show that after s successful purification 
steps the fidelity will be 

F s <l-(2/3) s (l-F )e . (19) 

In order to calculate Pq(s) we note that to obtain one pair of fidelity F s one uses up at least 2 S_1 pairs, and performs 
in each of them at least one operation. Therefore, Pq(s) > 2 S ~ 1 U, where now U denotes the joint cost per used pair. 
Upon substituting these expressions in ( |l8| ) we see that unless U is much smaller than the other costs or Fq ~ 1, there 
is no gain at all by using entangled states. 
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B. Purification Scheme 2 



We will evaluate now the cost associated with pre-computation when the nodes are networked via photonic channels 
||. As in the previous scheme, the maximum fidelity resulting in this case, F s , approaches 1 with the number s of 
purification steps exponentially fast 

F s ~ 1 - a s (l - F ), (20) 

where a < 1 is a constant. In this scheme the cost of establishing one entangled pair Po(s) — bs is proportional to 
the number of steps. Substituting these quantities in Eq. ( |l8| ) one can see the existence of an interval [n m , n , n max ] in 
which the use of maximally entangled states is more convenient with respect to uncorrelated ones. 

This purification protocol is more efficient than the previous one since Pq does not scale exponentially with s. The 
reason for this scaling law in the purification scheme 1 is that one discards one pair at each purification step, which 
makes it very ineffective. Therefore, it is more desirable to use a purification scheme giving a finite yield, as the 
hashing or breeding methods |J. Moreover, the purification protocols used here are just based on the ones developed 
for two qubits, which lead to the exponential dependence of F on n. By using ideas similar to the ones developed for 
quantum repeaters []l0|| it may be possible to improve this exponential dependence. 



VI. NOISE IN THE COMPUTATIONAL PROCESS 



Let us now assume that at each node the computation itself is not error-free but dephasing-type decohcrcncc is 
present at a rate <?, namely a random phase is introduced in front of the component | 1) of the qubit with probability 
e~ gt at time t (notice that if one considers a quantum optical implementation, the results that we will present in the 
following hold also in the presence of spontaneous emission). When measured in the computational basis, the bit 
value will now be obtained with probabilities 

p 1 = I[l-cos^-^)e-^] (21) 
when dealing with independent processors and 

P2 = - [1 - cos(n0 - ncj> 2 )e- n ^] + (22) 

in the case of maximally entangled nodes of the form (|J) . In the above equations t c is the time required to perform 
computation C and will be regarded as a fixed parameter in the comparison of the two schemes. To achieve resolution 
e the computation C must be performed 



ne 2 sin 2 (0 — (f>\) 



R 1 = -^ c ," 2 7, (23) 



times with independent processors and 



1 e 2^c f 2 cos 2(^_„ fe ) 

H z- n 2 e 2 F 2 sin 2 (n^-n0 2 ) 1 ' 

times when a maximally entangled input is distributed among the n nodes. The relative cost of both procedures is 
no longer <f> independent. Following the argument presented in the previous Section, we select the controllable phases 
<j)i^2 in both procedures in such a way that the measured phase approaches 7r/2 when dealing with uncorrelated inputs 
and 7r/2n when one uses entangled states. Therefore, we can write 

Rx = -^e 2 ^ (25) 



ne 2 



while 



The cost ratio is then given by 



i?2 = ^?2 e2 " 9tc - (26) 
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|=e-(-)(|) fl q . (27) 

As it can be seen from Figure 3, the effect of dephasing type decoherence is negligible in the limit gt c « 1/n. The 
net effect of decoherence during the computational process is a further reduction of the domain where the use of 
entanglement is cost efficient, being n max (g ^ 0) < n max (g = 0). 



VII. CONCLUSIONS 



In summary, we have introduced the notion of a generic cost of physical operations in distributed quantum com- 
putation. This parameter allows to quantify the efficiency of a quantum computation which can be run separately 
on different quantum processors belonging to a quantum network. Previous work [p]JTl|] has shown that the use of 
entangled states could be advantageous for certain computations. However, it was not obvious that the cost of gener- 
ating entanglement or the inclusion of noise during the computational process, might not nullify their advantage. We 
have shown that under certain circumstances a quantum network supplied with a maximally entangled input yields a 
smaller global cost than the one required when dealing with n independent inputs, provided that n belongs to a given 
interval n £ {n m i n ,n max }. We have illustrated this for the case of phase estimation. The values of n m in and n max 
crucially depend on the purification protocols used to establish the n- processor entangled states, as well as on the 
presence of decoherence processes during the computation. 
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FIG. 1. Experimental set-up for performing a distributed quantum computation between three nodes sharing an entangled 
state. 



FIG. 2. Cost ratio C2/C1 as a function of the number of nodes within the quantum network for two different sets of values 
of the parameters (X, Y, Z) = (100, 10, 1), (1000, 10, 1). Above a threshold n t h ~ X/Y , the global cost of the computation is 
smaller when using a maximally entangled input. 
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FIG. 3. Cost ratio as a function of the number of nodes for a phase estimation problem. The dotted line corresponds to the 
case when dephasing type decoherence takes place during the computational process. The net effect is that the domain where 
the use of entanglement is cost efficient gets shrinked with respect to the case of error free computation, represented by the 
solid line. 
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